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1.  Introduction 


In  the  competing  risks  model,  a  unit  is  exposed  to  several  risks  at  the  same 
time,  but  it  is  assumed  that  the  eventual  failure  of  the  unit  is  due  to  only  one  of 
these  risks  which  is  called  a  ‘cause  of  failure.’  Let  a  unit  be  exposed  to  two  risks  and 
the  notional  (or  latent)  lifetimes  of  the  unit  under  these  two  risks  be  denoted  by 
A'  and  Y,  respectively.  In  general,  A'’  and  1'  are  dependent.  Also,  being  lifetimes, 
they  are  nonnegative.  We  only  observe  {T,6)  where  T  =  min(A'’,y')  is  the  time  of 
failure  and  6  =  2  —  /(A  <  A)  is  the  cause  of  failure.  Here  /(.4)  is  the  indicator 
function  of  the  event  A. 

On  the  basis  of  the  competing  risks  data  it  is  often  useful  to  distinguish  between 
the  following  alternatives;  (i)  the  forces  of  the  two  risks  are  equal,  and  (ii)  the  force 
of  one  risk  is  greater  than  that  of  the  other,  within  the  environment  in  which  the 
two  risks  are  acting  simultaneously.  Such  comparisons  can  be  made  in  terms  of 
cumulative  incidence  functions  and  cause-specific  hazard  rates,  defined  as  follows. 
The  cumulative  incidence  function  corresponding  to  cause  j  is 

Fj{t)  =  P{T  <t,6  =  j], 

and  the  cause- specific  hazard  rate  (CSHR)  for  cause  j  is 


9}{i)  =  fj{t)/ST{t), 

where  the  Fj  are  assumed  to  have  subdensities  fj{t),  and  5j’(t)  =  P{T  >  t]  = 
1  —  Fi{t)  —  F2it)  is  the  survival  function  of  T.  In  the  case  when  A'  and  Y  are 
independent,  g\  and  g2  reduce  to  the  hazard  rates  corresponding  to  the  marginal 
distributions  of  A'  and  Y.  Prentice  et  al.  (1978)  show  that  in  general  only  proba¬ 
bilities  expressible  as  functions  of  and  g2  may  be  estimated  from  the  observable 
data  (T,  6). 

We  mention  two  practical  examples  in  which  it  is  important  to  be  able  to  make 
comparisons  between  cause-specific  hazard  rates  (as  well  as  cumulative  incidence 
functions).  The  first  example  arises  in  reliability  testing.  Suppose  that  either 
of  two  components  in  a  series  system  can  be  replaced  to  improve  overall  system 
reliability.  One  would  replace  the  second  component  in  preference  to  the  first  if 
<  gi  ov  Fi  <  F2.  The  second  example,  this  one  biomedical,  comes  from  a  paper 
of  Benichou  and  Gail  (1990).  It  concerns  time  to  cancer  recurrence  in  patients 
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with  surgically  resected  cancer.  In  deciding  whether  to  give  a  toxic  therapy  in  the 
hope  of  preventing  cancer  recurrence  it  is  appropriate  to  compare  the  cumulative 
incidence  function  for  cancer  recurrence  (Benichou  and  Gail  call  this  the  absolute 
risk  of  recurrence),  Fi,  and  the  cumulative  incidence  function  for  other  risks  (e.g. 
other  causes  of  death),  F2.  A  physician  would  be  reluctant  to  recommend  a  toxic 
cancer  treatment  in  an  elderly  patient  in  whom  F]  <  F2  or  <52-  Benichou  and 
Gail  go  on  to  discuss  the  importance  of  the  concept  of  absolute  risk  in  evaluating 
public  health  measures  to  prevent  disease.  Gray  (19SS)  also  draws  attention  to  the 
usefulness  of  comparing  cause-specific  incidence  for  different  types  of  failure. 

In  this  article  we  propose  some  methods  for  comparing  CSHR’s  and  cumulative 
incidence  functions.  Our  tests  are  less  subjective  than  graphical  procedures  based 
on  inspections  of  estimates  of  the  CSHR’s  themselves.  We  are  interested  in  testing 
the  hypotheses 

Hi  :  gi{t)  <  g2{t),  t  >  0 

H2  :  Flit)  <  F2{t),  t>0 

with  strict  inequalities  for  some  t.  These  hypotheses  represent  different  ways  of 
saying  that  risk  Y  is  “more  serious”  than  risk  X.  Clearly  Hi  implies  H2.  We 
regard  Hi  and  H2  as  alternatives  to  the  null  hypothesis  of  equal  risks: 

Ho  :  9iit)  =  g2{i),  t  >  0. 

Note  that  there  is  often  no  reason  to  expect  a  vriori  that  the  cause  specific  risks  gi 
and  g2  are  equal  (except,  say,  when  gi  and  g2  represent  two  identical  components  in 
a  series  system),  but  this  is  the  natural  choice  of  null  hypothesis  against  which  the 
ordered  alternatives  Hi  and  H2  should  be  tested.  A  similar  choice  of  null  hypothesis 
is  made  in  the  two-sample  survival  analysis  problem  when  testing  whether  survival 
in  one  group  is  better  than  survival  in  another  group,  cf.  Pepe  and  Fle  ming  (19S9). 

In  the  case  that  A'  and  Y  are  independent,  Bagai,  Deshpande  an  a  Kochar  (19S9 
a,b)  proposed  distribution-free  tests  for  testing  the  equality  of  two  competing  risks 
against  stochastic  ordering  and  failure  rate  ordering.  Sen  (1979)  proposed  nonpara- 
metric  tests  with  maximum  asymptotic  relative  efficiency  for  interchangeability  of 
the  competing  risks  against  alternatives  specified  in  teims  of  P[6  =  j\T  =  t].  Aras 
and  Deshpande  (19S9)  p^-oposed  locally  most  powerful  rank  tests  for  testing  Ho 
against  various  parametric  alternatives  expressed  in  terms  of  Fi  and  F2. 
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The  reiiuiiiicler  of  the  paper  is  orgaiiized  as  follows.  la  Section  2  we  consider 
the  problem  of  testing  the  null  hypothesis  of  equality  of  CSHR's  against  the  alter¬ 
natives  specified  by  Hi  and  H2-  Our  tests  are  of  Kolmogorov-Smimov  type  and  axe 
distribution-free.  Formulae  for  the  exact  null  distributions  of  the  test  statistics  are 
provided.  The  asymptotic  distributions  are  also  derived.  In  Section  3  we  consider 
an  extension  of  our  tests  to  deal  with  right-censored  data.  The  results  of  a  simu¬ 
lation  study  and  an  example  are  discussed  in  Section  4.  In  Section  5  we  describe 
another  extension  of  our  tests  to  allow  comparisons  on  any  given  time  interval. 

2.  Testing  for  the  equality  of  CSHR’s  against  ordered  alternatives 

In  this  section  we  introduce  our  tests  for  H\  and  H2  in  the  uncensored  case.  The 
tests  are  based  on  the  competing  risk  data  {{T,,6,)\i  =  1.. . .  ,n}  for  n  independent 
and  identical  units. 


2.1  Testing  Hq  against  Hi.  A  distribution-free  test  of  Hq  vs.  Hi  can  be  constructed 
using  the  fact  that,  under  ff],  the  function  xp{t)  =  Fi{t)  —  F2{t)  is  nonincreasing 
in  t.  This  is  a  consequence  of  the  identity  Fj{t)  =  /J  gj{u)ST{u)  du  and  provides  a 
rationale  for  the  test  statistic 


sup  {0,d'n{s)  -  ll'n{t)}, 

0<a<t<oo 

where  d'nit)  =  Fi„{t)  -  F2„it).  Here  Fj„{i}  =  n~^  I{6,  =  j.  T,  <  t}  is  the 

empirical  estimator  of  Fj.  Positive  values  of  provide  evidence  that  92  is  larger 
than  gi  in  some  interval.  Note  that 


where 


Dtu  = 


1 

max  — 
0<i<2<n  n 


(=t  +  l 


1 

=  max  — 
0<i<2<n  n 

1 

=  max  — 

0<i<j<n  11 


J  1  if  6  corresponding  to  T(,)  (the  ith  ordered  T,)  is  1 
[  0  otherwise. 


3 


Ut  =  1  -  211’  ,  /  =  i,  2, . . . ,  Zo  =  0  and  Z/-  =  ;/]+...  +  7/^.,  h  =  1,2, - n. 

To  obtain  the  exact  null  distribution  of  we  argue  cls  follows.  Kochar  and 
Proschan  (1991)  proved  that  T  and  8  are  independent  under  Hq.  Consequently, 
under  i/o,  are  i.i.d.  Bernoulli  random  variables  with  P(1T,  =  0)  = 

P(1P-  =  1)  =  i  and  7?i,...,77„  axe  i.i.d.  with  P(77,  =  —  1)  =  P(r7,  =  1)  =  i.  It 
follows,  using  Simons  (19S3),  that 


<  t}  =  P{—t  <  min  Zj  <  max  Zj  <  i  +  1},  (2.2) 

0<j<n  0<j<n 

cf.  Alj',  Gombay  and  Kochar  (1992,  proof  of  Theorem  3.1  (i)).  By  applying  Csaki 
(19S6,  (2.6)),  this  yields  an  e.xact  formula  for  the  null  distribution  function  of 
The  asymptotic  null  distribution  can  be  obtained  from  (2.1)  and  the  discussion 
preceding  (2.2)  by  using  the  invariance  principle  for  partial  sums  (see.  for  example. 
Chapter  2  of  Csdrgo  and  Revesz  (19S1)).  We  summarize  this  discussion  in  the 
following  theorem. 


Theorem  2.1.  Under  Ho 


j=0 


(2.3) 


sm 


jt: 


2t  +  l 


for  t  =  l,...,7r  +  l,  and 


\/nP+„-^  sup  |IK(x)l 

0<I<1 

where  {ir(t),t  >  0}  is  a  standard  Brownian  motion.  Consequently,  for  c  >  0 

4 

k=0 


P{v^P+„  <c}  ^  ^J2^j-^exp{-7:^{2k  +  lf/Sc^}.  (2.4) 


The  exact  formula  (2.3)  can  be  used  to  generate  a  table  of  critical  values,  see 
Aly,  Gombay  and  Kochar  (1992).  Using  (2.4)  the  asymptotic  0.90,  0.95  and  0.99 
quantiles  of  \/nD'l^  are  found  to  be  1.96,  2.241  and  2. 807,  respectively. 


2.2  Testing  Hq  against  H2.  Since  under  H2-,  4'it)  is  nonpositive,  a  natural  test 
statistic  for  testing  Ho  against  H2  is  given  by 

=  sup  {0, -Vn(0}- 

0<f  <oc 
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We  reject  H{j  for  large  values  of  Note  that  D2,,  can  be  written  as  maxo<;<n  Zj 
where  Zo,  Zi,. . .  ,Z„  are  as  in  (2.1).  Thus,  by  Lemma  4.8.1  of  Renyi  (1970), 


=  =  I- =  0,1,2 . 
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under  Hq.  This  gives  the  exact  null  distribution  of  The  asymptotic  null 

distribution  is  obtained  using  the  invariance  principle  for  partial  sums:  under  Hq 


>x]  ^  P{  sup  W'(0  >  x)  =  2(1  -  4>(x)),  X  >  0, 

0<t<l 

where  ^  is  the  standard  normal  distribution  function. 


3.  Censored  data 


In  this  section  we  consider  an  extension  of  our  tests  to  allow  for  the  possibility 
of  right-censoring.  The  underlying  censoring  mechanism  will  be  represented  by  a 
censoring  time  C  which  is  assumed  independent  of  the  latent  failure  times  A'  and 
Y.  Denote  the  survival  function  of  C  by  Sc  and  assume  that  Sc(t)  >  0  for  all  t. 

Under  right-censoring  we  observe  n  iid  copies,  (Ti,Si),  i  =  1 _ ,n.  of  T  = 

min(r,  C)  and  S  —  SI(T  <  C).  Our  approach  is  to  look  for  a  suitable  modification 
of  the  function  i/j.  Recall  that  t/>(t)  =  J*  ST{u  —  ){gi{u)—g2{u))  du.  In  order  to  obtain 
distribution-free  tests  of  Hi  and  H2  in  the  censored  case,  look  at  the  function 

4>{t)=  j  ST{u-)Sc{u-y''^{g\{u)  -  g2[u))du. 

Jo 

The  factor  —  turns  out  to  be  precisely  what  is  needed  to  compensate  for 
the  censoring.  Wb  have  that  nonincreasing  under  Hi  and  identically  zero  under 
Ho,  so  to  test  Hq  against  Hi  it  is  natural  to  use  the  test  statistic 

=  sup  {Q,0n{s)  - 

0<  j<t<oo 

where  4>„  is  a  suitable  estimator  of  6.  Similarly, 

Z>4n  =  sup  {0, -(?„()■)}, 

0<«<oo 
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can  be  used  to  test  Ho  against  Ho-  An  obvious  choice  of  is 

Mt)=  /  ST{u-)Sc{u-y^^diA,-A2){u), 

Jq 

where  St  and  Sc  are  the  product-limit  estimators  of  St  and  Sc,  and  Aj  is  the 
Aalen  estimator  of  the  cumulative  CSHR  function  Aj{t)  =  gj(u)  du: 

Aj(t)=  IiS,=j)lR, 

where  R,  =  ^{k  :  Ti^  >  T,}  is  the  size  of  the  risk  set  at  time  T,  — . 

We  note  that  Aj  is  a  special  case  of  an  estimator  discussed  by  Aalen  and 
Johansen  (1978)  in  connection  with  estimation  of  the  transition  probabilities  of  a 
non-homogeneous  Markov  chain  with  finitely  many  states.  Indeed,  we  are  dealing 
with  a  three-state  non-homogeneous  Markov  chain  having  two  absorbing  states 
corresponding  to  the  two  types  of  failure. 

The  following  result,  proved  in  the  Appendix,  shows  that  and  are 
asymptotically  distribution-free  and  have  the  same  limiting  distributions  they  do  in 
the  uncensored  case. 

Theorem  3.1.  Under  Ho 

sup  |Ti’’(x)|  and  sup  lT(x). 

0<I<1  0<i<l 

Our  approach  easily  extends  to  the  case  of  multiple  (rather  than  just  two) 
competing  risks  in  which  any  two  of  the  cause-specific  risks  are  to  be  compared. 
No  structure  needs  to  be  imposed  on  the  dependency  between  the  multiple  risks, 
although  the  corresponding  latent  failure  times  need  to  be  independent  of  the  cen¬ 
soring,  as  before.  Let  T  be  the  minimum  of  a  finite  collection  of  latent  failure  times 
which  include  A'  and  1',  and  let  ^  denote  the  corresponding  cause  of  failure.  Exten¬ 
sions  of  and  tnat  preserve  the  above  asymptotic  distributions  are  obtained 
by  using  (i>n{i)/y/p^  in  place  of  dn(0,  where 

rOO 

Pn  =  /  St{u-)  d{Ai  +  A2){u) 

Jo 

is  a  consistent  estimator  of  P [6  =  1  or  2],  see  the  Appendix. 
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4.  Simulation  x'esults  and  an  example 


Our  test  procedures  are  consistent  against  their  respective  alternatives,  H\  or 
H^-  However,  we  would  like  to  know  whether  they  are  powerful  enough  for  practical 
applications.  For  that  purpose,  we  carried  out  a  simulation  stud}',  and  the  results 
show  that  our  tests  are  readily  able  to  detect  these  ordered  alternatives. 

For  the  distribution  of  (.Y,  K)  we  used  Block  and  Basu's  (1974)  absolutely 
continuous  bivariate  exponential  (ACBVE)  distribution  having  density 

(  Ai A(A;+Ao)  -Aii-(A;  +  Ao)v  -y  c'  t, 

f(^  ^  J  A,+A,  ^  nu,<y 

I  ■^?-^(-^i  +  -^p)^-A?v-(Ai+An)i  X  <  y 

where  (Aq,  Aj  ,  A2 )  are  parameters  and  A  =  Aq  +  Ai  +  A2.  The  CSHR's  are  given  by 


9j{i)  = 


AjA 


-^1  +  ^2 

so  Hi  holds  if  and  only  if  Aj  <  A2.  Under  this  model  Hi  and  H2  are  equivalent. 
The  parameter  Aq  controls  the  degree  of  dependence  between  X  and  Y;  they  are 
independent  if  and  only  if  Aq  =  0.  We  set  Aj  =  1  and  considered  various  higher 
values  of  A2  corresponding  to  larger  and  larger  departures  from  Ho-  The  censoring 
was  taken  to  be  exponential  with  parameter  values  1  and  3,  corresponding  to  “light’’ 
and  “heavy”  censoring  (about  25%  and  50%  censored,  respectively).  For  the  sake  of 
comparison  we  included  results  for  the  uncensored  case  as  well.  We  used  asymptotic 
critical  levels  of  5%. 

Inspection  of  Table  1  shows  that  use  of  the  asymptotic  critical  levels  gives 
somewhat  conservative  tests,  and  this  effect  increases  as  the  censoring  becomes 
more  severe.  The  test  based  on  appears  to  be  more  conservative  than  the 
one  based  on  However,  the  tests  become  less  conservative  as  the  sample 

size  increases  (in  fact  we  have  found  that  the  levels  of  the  tests  are  close  to  their 
nominal  5%  values  for  sample  sizes  over  500,  even  under  heavy  censoring).  There 
is  no  apparent  adverse  effect  on  the  levels  or  the  power  due  to  lack  of  independence 
of  A'  and  Y.  (Pearson’s  correlation  between  A'  and  Y  is  about  .15  for  the  table 
entries  corresponding  to  Aq  =  1.) 

As  an  application  we  have  analyzed  a  set  of  mortality  data  given  in  Hoel  (1972). 
These  data  were  obtained  from  a  laboratory  experiment  on  99  RMF  strain  male  mice 
which  had  received  a  radiation  dose  of  300  rads  at  5-6  weeks  of  age  and  were  kept 


in  a  couveatioual  laboratory  environment.  The  cause  of  death  was  classified  int  j 
thymic  lymphoma,  reticulum  cell  sarcoma,  and  other  causes.  For  us,  “other  causes” 
represents  censoring  (39%  were  censored),  and  the  two  types  of  cancer  mortality 
are  taken  to  be  the  two  causes  of  failure  that  we  wish  to  compare,  i.e.  gi  and  §2 
are  the  CSHR’s  from  lymphoma  and  sarcoma  respectively.  Our  analysis  depends 
on  the  assumption  that  the  two  diseases  are  lethal  and  independent  of  other  causes 
of  death,  but  we  do  not  need  to  assume  that  they  are  independent  of  one  another. 

[Insert  Figures  1  and  2  about  here] 

Figure  1.  Aalen  estimates  of  cumulative  CSHR’s  for  Ij'mphoma  (dashed  line)  and 
sarcoma  (solid  line). 

Figure  2.  Plot  of  \/no„(t)  (solid  line)  and  the  asymptotic  5%  critical  levels  for 
(dashed  lines). 

We  obtained  =  4. SI,  %vhich  gives  a  P-value  of  less  than  .01  for  testing  Hq 

against  Hi.  Also,  =  2.77,  which  gives  a  P-value  of  .0056  for  testing  H{j  against 

H2-  When  the  roles  of  lymphoma  and  sarcoma  were  reversed,  we  obtained  = 
=  2.03,  so  the  P-values  for  the  two  tests  are  close  to  .1  and  .05  respectively. 
Our  conclusion  is  that  the  two  cause-specific  hazard  rates  are  unequal.  Note 
that  we  cannot  conclude  that  the  CSHR  for  sarcoma  is  uniformly  larger  than  the 
eSHR  for  lymphoma;  the  large  value  =  4.81  only  indicates  that  the  sarcoma 
CSHR  is  larger  than  the  lymphoma  CSHR  in  some  age  interval.  Indeed,  inspection 
of  a  plot  of  the  two  cumulative  CSHR  estimates  (Figure  li  suggests  that  up  to  500 
days  there  is  moderate  risk  of  lymphoma,  yet  negligible  risk  of  sarcoma.  After  500 
days  the  situation  reverses:  there  is  negligible  risk  of  lymphoma  but  high  risk  of 
sarcoma,  and  it  is  this  large  difference  that  the  test  statistic  is  picking  up. 
This  is  also  reflected  in  the  plot  of  \/n0„(f)  in  Figure  2.  Such  plots  are  useful  in 
avoiding  misinterpretation  of  the  test  statistics.  Plots  of  estimates  of  the  CSHR’s 
themselves  are  also  useful;  these  can  be  made  by  finding  smoothed  derivatives  of  the 
cumulativ'e  CSHR  estimates,  see  Ramlau-Hansen  (19S3),  and  are  somewhat  easier 
to  interpret  than  plots  of  the  cumulative  CSHR  estimates.  However,  our  tests  offer 
a  less  subjective  comparison  of  CSHR’s  than  can  be  made  fren  a  simple  visual 
inspection  of  such  plots. 
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5.  Comparing  CSHR’s  in  [^1,^2) 

It  is  often  useful  to  compare  CSHR’s  (or  cumulative  incidence  functions)  in  a 
given  time  interval  [<1,^2)^  rather  than  at  all  times.  For  instance,  an  examination 
of  plots  of  the  cumiilative  CSHR  estimates  for  Hoel’s  data  strongly  suggests  that 
the  CSHR  for  sarcoma  is  much  larger  than  the  CSHR  for  lymphoma  after  500  days. 
In  the  second  example  discussed  in  the  Introduction,  Benichou  and  Gail  (1990,  p. 
820)  are  interested  in  comparing  the  CSHR  for  cancer  recurrence  with  the  CSHR 
for  other  risks  at  times  between  one  and  five  years  following  surgical  tT-eatment. 

It  is  straightforward  to  generalize  our  tests  in  Section  3  to  deal  with  such  cases. 
We  want  a  test  of 

Ho  :  giit)  =  g2{t),  <t  <t2 

against  the  alternative 

Hi  ;  gi{t)  <  g2{t),  h  <  t  <  t2 

with  strict  inequality  for  some  t  €  [^1,^2)-  We  replace  <f)  by  the  t’unction 

<?*(t)  =  (5r(ti)-5T(t2))"''''  /  5r(w-)5c(u-)'/'(5i(u)-^2(t^))rfu. 

Jtx 

Clearly  H*  is  equivalent  to  d*  nonincreasing  on  [ti,t2)-  As  before,  we  suggest  the 
test  statistic 

sup  {0,<^;(5) 

ti  <S<<<<2 

where  d*,(t)  is  obtained  by  substituting  St  etc.  into  d*.  It  can  be  shown  by  routine 
modifications  of  the  proof  of  Theorem  3.1  that  converges  in  distribution  to 

siiPo<x<i  under  Hq. 

When  this  test  was  applied  to  Hoel’s  data,  we  obtained  the  highly  significant 
values  of  =  5.56  (resp.  3.69)  when  testing  whether  the  CSHR  for  sarcoma  is 
larger  (resp.  smaller)  than  the  CSHR  for  lymphoma  after  (resp.  before)  500  days. 
This  confirms  our  earlier  conjectuT-es  arising  from  examination  of  Figures  1  and  2. 

APPENDIX 

Proof  of  Theorem  3.1.  Suppose  we  can  show  that 

VV^d'n^WiFri-)).  (A.l) 
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Then,  the  second  part  of  the  theorem  is  clear.  Using  the  continuous  mapping 
theorem, 


sup  sup  {W''(i^r(s))-TU(Fr(0)) 

0<s<<<oc  0<s<<<oo 

=  sup  {TU(n)  —  T'U(t’)} 

0<u<r<l 

=  sup  {W{u)  —  inf  TU(u)}. 

0<«<1  0<r<u 

The  statement  of  the  first  part  of  the  theorem  now  follows  by 
W{u)-  inf  lU(u)=  sup  ir(u)  - 

0<»’<“  0<f<ii 


and  the  following  well  known  result  of  Levy  (194S): 


sup  T'l'fe)  —  lU(a)=]IU(u)|, 

U  <  i'  <  ti 


see  Chung  and  Williams  (19S3).  It  remains  to  prove  (A.l),  for  which  we  use  the 
counting  process  approach  developed  by  Aalen  (1978).  Note  that  we  can  write  Aj 
in  the  form 


dNjiu) 
F(u)  ’ 


where  1/0  =  0, 

r,(u)  =  Hf.  >  u),  N„(u)  =  I(f,  <  u,h  =  j), 


for  j  —  1,2,  and  the  summations  are  over  z  =  1, . . . ,  n.  Let 

Jo 

Then  Mij,i  =  1, ...  ,n  are  orthogonal  martingales  under  the  natural  filtration  gen¬ 
erated  by  the  above  processes.  Let  A/j  =  ^  M,y  .  The  predictable  variation  process 
of  I\Ij  is  y(zi)  dAj(u).  By  P{X  =  F)  =  0,  the  counting  processes  A'l  and 
almost  surely  have  no  simultaneous  jumps,  so  A/j  and  M2  are  orthogonal  martin¬ 
gales  (this  is  a  standard  resuT  from  counting  process  theory).  Thus,  the  predictable 
variation  process  of  Mi  —  M2  is  Jj  Y{u)  dAo{u),  where  Aq  =  Ai  -|-  Ao.  Under  Hq 


<Pn{t) 


I 


'  5r(u-)5c(u-)»/- 
Yin) 


d{Mi  -Yl2){u). 
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Since  S''/’(u  — )  and  Sc(a  — )  are  left  continuous  and  adapted,  they  are  predictable, 
so  \/n<^„  is  a  martingale  with  predictable  variation  process 


5c(a-) 

)  n^yn 


dAo{u). 


By  the  Glivenko-Cantelli  theorem,  Y{u)fn  converges  uniformly  in  u  to  P{T  >  u)  = 
ST{'u  —  )Sciu  —  )  almost  surely.  Hence,  by  the  uniform  consistency  of  the  product- 
limit  estimator  on  [0,t],  the  above  variation  process  converges  in  probability  to 
fo  ST(u-)dAo  (u)  =  FT(t).  Here  we  have  used  the  fact  that  the  cumulative  hazard 
function  of  T  is  Aq;  see  Prentice  et  ah  (1978).  The  appropriate  Lindeberg  condi¬ 
tion  is  easily  checked.  (A.l)  follows  by  Rebolledo’s  (1980)  martingale  convergence 
theorem.  □ 


We  conclude  by  indicating  how  to  extend  the  above  proof  to  deal  with  multiple 
competing  risks.  In  this  setting  the  predictable  variation  process  of  \/n<f>n  converges 
in  probability  to  Fj  -f-  Fz.  Since  pn  is  consistent  for  P[S  =  1  or  2],  it  follows  that 

V  Pn 

where  F\2  is  the  conditional  distribution  function  of  min(A*,  V')  given  that  6  =  1  or 
2.  This  extends  (A.l).  The  remaining  steps  of  the  proof  are  identical. 
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Table  1.  Observed  levels  and  powers  of  tests  for  equality  of  CSHR's  based  on 
(resp.  D^n)  at  an  asymptotic  level  of  5%.  The  underlying  distribution  of  (A',}'')  is 
Block  and  Basu’s  (1974)  ACBVE  with  Ai  =  1. 

(b)  Uncensored 


X2 

n  =  50 

n  =  100 

0 

II 

0 

Ao  —  1 

0 

II 

0 

Ao  —  1 

1.0 

3.86  (4.90) 

3.86  (4.90) 

3.68  (4.44) 

3.69  (4.44) 

1.5 

32.37  (39.46) 

32.38  (39.46) 

54.41  (61.05) 

54.43  (61.05) 

2.0 

67.46  (74.95) 

67.46  (74.95) 

92.59  (95.11) 

92.59  (95.11) 

2.5 

87.66  (91.96) 

87.66  (91.96) 

99.4  (99.78) 

99.40  (99.78) 

(b)  Lightly  censored  (18%-33%) 


A2 

n  =  50 

n  =  100 

Ao  —  0 

Ao  —  1 

0 

II 

0 

Ao  —  1 

1.0 

2.89  (3.64) 

2.98  (3.87) 

3.45  (4.16) 

3.37  (4.06) 

1.5 

21.95  (27.64) 

24.19  (30.00) 

41.32  (47.97) 

44.86  (51.22) 

2.0 

51.95  (60.52) 

55.40  (63.64) 

83.31  (87.64) 

85.57  (89.76) 

2.5 

76.35  (82.91) 

78.49  (84.80) 

97.47  (98.57) 

97.97  (98.75) 

(c)  Heavily  censored  (40%-60%) 


A2 

n  =  50 

n  =  100 

0 

II 

0 

Ao  —  1 

Ao  =  0 

Ao  —  1 

1.0 

1.49  (2.29) 

2.16  (2.82) 

1.80  (2.61) 

2.79  (3.64) 

1.5 

11.09  (16.02) 

14.76  (19.79) 

22.88  (29.12) 

29.24  (35.85) 

2.0 

30.38  (39.76) 

37.26  (46.75) 

60.59  (68.79) 

69.25  (76.49) 

2.5 

53.11  (63.73) 

60.73  (70.27) 

86.49  (91.57) 

91.01  (94.72) 

Note;  The  data  were  created  using  the  uniform  random  number  generator  of  Marsaglia, 
Zaman  and  Tsang  (1990)  and  an  algorithm  of  Friday  and  Patil  (1977,  Corollary  3.3). 
10000  samples  were  used  to  obtain  each  entry  in  the  table. 
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